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ABSTRACT 

Context. Over the past few years several studies have provided estimates of the SFR (star-formation rate) or the total infrared luminosity from just 

one infrared band. However these relations are generally derived for entire galaxies, which are known to contain a large scale diffuse emission that 

is not necessarily related to the latest star-formation episode. 

Aims. We provide new relations to estimate the SFR from resolved star-forming regions at 100 fim and 160 yum. 

Methods. We select individual star-forming regions in the nearby (840 kpc) galaxy M33. We estimate the SFR combining the emission in Ha 

and at 24 /um to calibrate the emission at 100 yum and 160 /im as SFR estimators, as mapped with FACS/Herschel. The data are obtained in the 

framework of the HERM33ES open time key project. 

Results. There is less emission in the HII regions at 160 /jm than at 100 fim. Over a dynamic range of almost 2 dex in IIsfr we find that the 100 /jm 

emission is a nearly linear estimator of the SFR, whereas that at 160 yum is slightly superlinear 

Conclusions. The behaviour of individual star-forming regions is surprisingly similar to that of entire galaxies. At high 1.sfr, star formation 

drives the dust temperature, whereas uncertainties and variations in radiation-transfer and dust-heated processes dominate at low 1.sfr. Detailed 

modelling of both galaxies and individual star forming regions will be needed to interpret similarities and differences between the two and assess 

the fraction of diffuse emission in galaxies. 
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1. Introduction 

Star formation is one of the main drivers of galaxy formation and 
evolution and as such the accuracy of the SFR (star-formation 
rate) determination is of great importance for deriving the cos- 
mic history of galaxies. Along with the UV and the Ha, the to- 
tal infrared luminosity is widely used to estimate the SFR. To 
properly quantify the infrared luminosity a good sampling of 
the infrared SED (spectral energy distribution) is needed (Dale 
& Helou 2002; Draine & Li 2007). Over the past few years, 
many authors have shown that the TIR (total infrared) luminos- 
ity, and by extension the SFR (Kennicutt 1998), can also be eval- 
uated from monochromatic emission measures (Takeuchi et al. 
2005; Calzetti et al. 2007; Rieke et al. 2009; Calzetti et al. 2010; 
Boquien et al. 2010, Y. Li et al. 2010, in prep.). However, most 
relations between the far-infrared luminosity and the SFR are 
established for entke galaxies. At shorter wavelengths, some re- 
lations have been derived from individual star-forming regions 
(Calzetti et al. 2005; Perez-Gonzalez et al. 2006; Calzetti et al. 
2007; Relano et al. 2007). As new deep surveys will become 
available at wavelengths where most of the energy is emitted, it 
is prudent to try to understand the relation between the integral 
emission of a galaxy and that of the individual star forming re- 
gions. The physical conditions, such as temperature, abundance 
and emissivity, of the infrared-emitting dust can vary widely in 
a galaxy, and so does the significant contribution from evolved 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



Stars (Lonsdale Persson & Helou 1987; Sauvage & Thuan 1992; 
Buat & Xu 1996). Therefore, any scaling relation established 
from the emission of entire galaxies may not be appropriate 
when applied to resolved star-forming regions in these same 
galaxies. 

The limited resolution of far-infrared instruments onboard 
IRAS (Infrared Astronomical Observatory), ISO (Infrared Space 
Observatory) or even Spitzer beyond 60 pm made the study of in- 
dividual star-forming regions within galaxies difficult. However, 
the recently launched Herschel Space Observatory (Pilbratt et al. 
2010) with its unprecedented resolution provides the first oppor- 
tunity to study the spatially resolved far-infrared dust emission 
in exquisite detail. Such a fine resolution is of the utmost impor- 
tance to study the emission of star-forming regions located in 
nearby galaxies in order to provide insights into the fundamental 
properties of the dust and to quantify the SFR. 

With an inclination of 56° (Regan & Vogel 1994) and a dis- 
tance of only 840 kpc (Freedman et al. 1991), M33 is one of 
the closest spiral galaxies. It has been imaged by Herschel in 
the context of the HERM33ES Open Time Key Project (Kramer 
et al. 2010), providing one of the finest views Herschel will ever 
provide of a spiral galaxy from 100 pm to 500 pm (Kramer et al. 
2010; Braine et al. 2010; Verley et al. 2010b). 

2. Observations and data reduction 

2.1. PACS 

The observations provided by PACS (Poglitsch et al. 2010) at 
100 pm and 160 pm are presented by Kramer et al. (2010) along 
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with a detailed description of the data processing pipeline. The 
observations were carried out on 2010-01-07 in parallel mode 
with a 20"/s scanning speed for a total of 6.3 hours, through 
a single scan and a perpendicular cross scan. The frames were 
first processed to level 1 with HIPE (Ott 2010), the drifts were 
corrected and were deglitched with the second-order deglitcher 
with a 6-cr threshold. The maps were produced with photpro- 
ject mapmaker' using a two-step masking technique to preserve 
the diffuse emission from being affected by the high-pass filter 
The total flux of the galaxy is consistent with the measures pro- 
vided by IRAS and Spitzer at 100 /im and 160 /vm (Kramer et al. 
2010). We present the two maps in Figure 1. The pixel size is 
3.2" at 100 ;um and 6.4" at 160 /i/m for a spatial resolution of 
6.7"x6.9" at 100 jum and 10.7"xl2.1"at 160 //m. The absolute 
calibration uncertainty is 5% at 100 //m and 10% at 160 yum. The 
total fluxes of M33 agree to within a few percent with those from 
ISO (Hippelein et al. 2003) and MIPS. In addition, for all radial 
averages the PACS 160 flux is within 20% of the MIPS 160 one. 

2.2. Ha and Spitzer MIPS 24 nm 

We used the Ha image presented by Hoopes & Walterbos (2000) 
that is commonly used in the literature (Tabatabaei et al. 2007; 
Gardan et al. 2007; Verley et al. 2007, 2009, 2010a). The MI 
contamination was corrected assuming [NIIl/Ho-^O.OS in the fil- 
ter bandpass. We also corrected the fluxes for Galactic fore- 
ground extinction using the Cardelli et al. (1989) law, assuming 
E(B-V)=0.042 from the NASA Extragalactic Database-. 

We used the 24 jum MIPS data presented by Tabatabaei et al. 
(2007). No further processing was performed on this image. 

2.3. Flux measurements 

All targeted HII regions flux densities were measured in polyg- 
onal apertures using IRAF's polyphot procedure. Each polygon 
was constructed manually from the PACS 160 /im image and 
tailored to avoid subtraction artefacts in Ho- and background 
sources in MIPS 24 yum images. Each source was selected to 
be as compact as possible, taking into account the blending at 
160 jum to avoid the mix of several star-forming regions of dif- 
ferent ages and properties. The background was calculated mea- 
suring the mode of the pixels distribution in an annulus around 
the aperture. Annulus pixels falling into the aperture of a source 
were automatically discarded. The inner radius of the annulus 
ranges from 20" to 70" by steps of 10", which were defined to 
be larger than the equivalent radius of the aperture: r = yfSjn, 
where S is the area of the aperture. A scale of 50 pc corresponds 
to an angular size of 11". The width of the annulus was set to 
12". Aperture correction was performed for 24 //m data^^ and for 
the PACS bands"*. As the apertures are not circular, we applied 
the method presented in Boquien et al. (2007) using the equiva- 
lent radius of the aperture. 

As a proxy for the SFR we applied the scaling presented by 
Calzetti et al. (2007): SFR = {L{Ha) + 0.03 IL (24)]x5.35 lO^^^ 
in M0 yr"', where L (Ha) is the Ha luminosity in W and, L (24) 



' PACS photometer - Prime and Parallel scan mode release note. 
V.1.2, 23 February 2010. 

^ The NASA/IPAC Extragalactic Database (NED) is operated by the 
Jet Propulsion Laboratory, California Institute of Technology, under 
contract with the National Aeronautics and Space Administration. 

' Following the formula provided in the MIPS instrument handbook. 

* Following the correction provided in PACS photometer - Prime and 
Parallel scan mode release note. V.1.2, 23 February 2010. 



is defined as vLy at 24 yum in W, assuming a Kroupa (2001) IMF 
(initial mass function) with a constant SFR over 100 Myr 

For an easier use of the SFR estimator we will provide in 
this article, we subsequently worked in E (luminosity suiface 
density) because it is distance-independent, in order to facili- 
tate a comparison with other galaxies. To do so, we divided the 
luminosity by the area of the polygon measured in kpc^. 



3. Results 

3.1. General characteristics of Hi! regions 

We selected a total of 179 HII regions from the 160 yum map. 
The physical equivalent radius of the extraction apertures ranges 
from 37 pc to 256 pc, with a median of 99 pc. Sioo ranges 
from 1.0 X 10^3 W kpc^^ to 1.0 x 10^^ W kpc^^ and Sieo from 
1.5 X 10^^^ W kpc-2 to 5.2 X lO^-* W kpc^^. The typical 1-cr 
uncertainties are 0.09, 0.06, 0.02, and 0.03 dex in Ha, 24 yum, 
100 fim, and 160 yum respectively. 'Lsfr (SFR density) ranges 
from 1.2 X 10"^ Mq kpc"^ yr^^ to 1.5 x 10"' Mq kpc^^ yr'. 

The fraction of the total flux enclosed in the 179 apertures 
compared to the total flux of M33 is 0.40, 0.43, 0.35, and 0.24 
in Ha, 24yum, 100 //m, and 160 //m assuming galaxy-integrated 
fluxes F(Ha) = 4.03 x 10"'^ W m^^, F(24) = 49.4 Jy (Verley 
et al. 2007), F(IOO) = 1288 Jy, and F(160) = 1944 Jy (Kramer 
et al. 2010). That so little of the total 160 yum flux associated 
with the selected HII regions suggests that the large-scale diff'use 
emission seen at this wavelength may not be directly related to 
the ongoing massive star-formation, which is consistent with the 
result of Hinz et al. (2004) for M33, but may be heated by non- 
ionising B and A stars as observed for instance by Israel et al. 
(1996) in another galaxy with a similar metallicity, NGC 6822. 



3.2. 100 fim and 160 fim as SFR estimators 

Being closer to the peak IR emission, the 100 //m promises 
to yield an accurate SFR indicator for HII regions and star- 
formation-dominated galaxies. In this respect, the high angular 
resolution Herschel data will be essential to establish the range 
of applicability and any limitation of this indicator. 

First of all Zioo and 2i6o are very well correlated with a 
Spearman correlation coefficient pioo-ieo - 0.92. This is ex- 
pected because both bands probe the grey body emission of big 
grains. It appears that Eioo and Ei6o are also well coiTelated with 
'^SFR, with a Spearman coiTelation coefficient psFR-\Qa - 0.85 
and psfR_ 160 = 0.77 . 

In Figure 2 we present the fits of Sioo and Sigo versus the es- 
timated "Lsfr- To estimate the relations between the SFR and the 
PACS emission we fitted a linear relation in log-log using an or- 
dinary least-square technique taking into account uncertainties 
on both axes. 

The best fits for M33 HII regions respectively correspond to 

logSsfR = (0.9938 ±0.0244)logSioo- (35.8520 + 0.8304X1) 
log^sFR = (1. 1333 ±0.0293) log Sieo- (40.4325 ±0.993 1X2) 

The relations have a dispersion around the best fit of 0.22 
dex and 0.25 dex. We notice that the 100 /im relation is nearly 
linear, whereas the 1 60 /im relation is slightly superlinear. This is 
similar to what Calzetti et al. (2010) found (Eq. 23), even though 
individual star-forming regions probe a much smaller range of 

'^SFR- 
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Fig. 1. PACS maps of M33 at 100 fim (left) and 160 jjm (right). To probe the full dynamic range of the image a logarithmic scale is 
used. Individual clumps are well resolved with minimal blending. Some diffuse emission is visible, particularly in the arms and the 
central region of the galaxy. The visible perpendicular stripes in each scanning direction are sampling artefacts. The red polygons 
represent the selected apertures. North is up. East is left 
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Fig. 2. 1,sFR as a function of Zjoo (left) and Zj^o (right). A black plus corresponds to individual HII regions in M33. The error bars 
correspond to the typical uncertainty for the lower and higher 50% percentile. For comparison additional integrated galaxies are also 
plotted: SINGS {Spitzer Infrared Nearby Galaxies Survey, Dale et al. 2007, red crosses), LVL (Local Volume Legacy Survey, Dale 
et al. 2009, green squares), starburst galaxies from Engelbracht et al. (2008) (blue triangles) and LIRGs from (Luminous Infrared 
Galaxies Calzetti et al. 2010, cyan circles). These galaxies belong to the high-metallicity bin of Calzetti et al. (2010), which 
encloses the metallicity of M33 (Magrini et al. 2010). The best fit of selected star-forming regions in M33, taking into account the 
uncertainties on both axes, is plotted with a solid black line. The dashed line represents the model described in Sect. 4. 1 assuming an 
extinction E(B-V)=0.25, set to reproduce the observed "Ena/'^iA- Finally the dotted line represents the best fit determined by Calzetti 
et al. (2010) for entire galaxies. 



4. Discussion 

4.1. Modeling 

To model the individual HII regions we used the Calzetti et al. 
(2007) model as a baseline. The ionising flux is determined using 
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Starburst99 (Leitherer et al. 1999) with a Kroupa (2001) IMF, 
an instantaneous burst and solar metallicity. The extinction is as- 
sumed to follow the Calzetti (2001) law. As the L{Ha)/L{24) 
ratio does not show a significant correlation with I.sfr for 
our selected regions in M33, we assumed a constant E(B- 
V)=0.25 mag, set to reproduce the observed mean l,Ha/'^24- For 
the dust emissivity we assumed the Draine & Li (2007) model 
prescriptions for Spitzer MIPS 160, pending updated dust emis- 
sivities for Herschel bands. The model is plotted in Figure 2. 

4.2. Comparison with entire galaxies 

Calzetti et al. (2010) showed that for entire galaxies the emission 
at 160/im correlates linearly with I.sfr- To compare this with the 
emission of entire galaxies at 100 yum and 160 fim, we plotted in 
Figure 2 several samples of relatively metal-rich galaxies from 
SINGS (Spitzer Infrared Nearby Galaxies Survey), LVL (Local 
Volume Legacy Survey), starburst galaxies from Engelbracht 
et al. (2008) and some LIRGS (Luminous Infrared Galaxies). 
Surprisingly, individual star-forming regions and entire galax- 
ies show a similar behaviour. Unfortunately the large scatter in 
both samples makes any assessment of the contamination by the 
large diffuse emission in entire galaxies very difficult. However, 
when they are compared with the relation derived for individual 
regions, most show a stronger 100 pm and 160 pm emission for 
a given 1.sfr, especially in the 33 < logE^ ^ 34.5 range where 
X is 100 or 160. Further comparisons to estimate the contami- 
nation by a diffuse far-infrared component will require detailed 
modelling of individual star-forming regions and entire galax- 
ies. 



4.3. Dust temperature 

The slight non-linearity in the Eigo - "Esfr relation hints at a 
higher dust temperature with increasing "Esfr, which would have 
different effects on the two PACS bands because the peak of the 
emission passes through the filter bandpasses as the tempera- 
ture increases. Indeed, as dust gets warmer, an increasing frac- 
tion will be emitted at shorter wavelengths. In Figure 3 we plot 
Si6o/2ioo, a proxy for the dust temperature of the warm compo- 
nent, versus the I.sfr- 

We see a clear trend with higher 1,sfr (p - -0.55) lead- 
ing to a higher dust temperature. The spanned range is compati- 
ble with the emissivity values published by Draine & Li (2007) 
for 0.5 < U < 100, U being the interstellar radiation field nor- 
malised to that of the solar neighbourhood. This means that an 
increasing fraction of the total dust emission is detected in the 
100 pm band compared to the 160 pm band. Interestingly we 
also notice that the trend in entire galaxies is very similar to the 
trend in individual HII regions in M33 despite the fact that indi- 
vidual regions should have little contamination from the diffuse 
large scale emission. One possible interpretation is that in both 
entire galaxies and individual HII regions, star-formation domi- 
nates at higher 1,5 fr and creates the trend whereas, at lower 'Esfr 
the trend is influenced by the uncertainties and variations of con- 
ditions in the radiation-transfer and dust-heating processes such 
as the opacity of the star-formation region, the dumpiness of the 
media, the relative locations of stars and dust clouds, etc. 



5. Conclusions 

We used the high-resolution Herschel 100 pm and 160 pm ob- 
servations of a nearby star-forming galaxy, M33. We combined 
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Fig. 3. Eiso/Sioo vs. 1,sfr (black plusses), for the 179 HII re- 
gions selected in M33. The error bars correspond to the typical 
uncertainty for the lower and higher 50% percentile. For compar- 
ison additional entire galaxies are plotted using the same colour 
scheme as in Fig. 2. 



Herschel PACS data with Spitzer MIPS 24 pm and ground- 
based Hff to provide new calibrations of the 100 pm and 160 pm 
to estimate 'Esfr from individual star-forming regions. For the 
selected star-forming regions in M33, the 100 pm luminosity is 
a linear SFR estimator over a factor 100 in surface brightness, 
whereas the 160 micron luminosity is slightly superlinear. It ap- 
pears that individual star forming regions exhibit a similar be- 
haviour as entire galaxies taken from the LVL, SINGS, starburst 
galaxies from Engelbracht et al. (2008) samples when estimating 
"LsFR from the 100 pm and 160 pm bands emission. In a similar 
fashion, the dust temperature - as measured by the ratio of the 
160 pm to 100 pm emission - increases as a function of I,sfr, 
suggesting that at high I,sfr the star formation drives the trend 
for both systems, while at lower 1,sfr uncertainties and vari- 
ations of conditions in the radiation-transfer and dust-heating 
processes contribute to the scatter. In other words, the fairly wide 
dust temperature distribution at low Esfr becomes increasingly 
biased towards higher temperatures at higher "Esfr in both HII 
regions and entire galaxies. 
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